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Abstract 

Motivated by the Kerr-CFT conjecture, we investigate perturbations of the near-horizon 
extreme Kerr spacetime. The Teukolsky equation for a massless field of arbitrary spin is solved. 
Solutions fall into two classes: normal modes and traveling waves. Imposing suitable (outgoing) 
boundary conditions, we find that there are no unstable modes. The explicit form of metric 
perturbations is obtained using the Hertz potential formalism, and compared with the Kerr- 
CFT boundary conditions. The energy and angular momentum associated with scalar field and 
gravitational normal modes are calculated. The energy is positive in all cases. The behaviour 
of second order perturbations is discussed. 
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1 Introduction 

Some time ago, Bardeen and Horowitz (BH) showed that one can take a near-horizon limit of the 
extreme Kerr geometry to obtain a spacetime similar to AdS2 x [1 . This near-horizon extreme 
Kerr (NHEK) geometry has an SL(2,R) x U{1) isometry group, where the U{1) is inherited from 
the axisymmetry of the Kerr solution and the SL{2, R) extends the Kerr time-translation symmetry. 
Recently, Guica, Hartman, Song and Strominger (GHSS) have conjectured that quantum gravity 
in the NHEK geometry with certain boundary conditions is equivalent to a chiral conformal field 
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theory (CFT) in 1+1 dimensions Using this, they gave a statistical calculation of the entropy 
of an extreme Kerr black hole. 

More precisely, GHSS showed that there exist boundary conditions on the asymptotic behaviour 
of the metric such that the asymptotic symmetry group is generated by time translations plus a 
single copy of the Virasoro algebra, the latter extending the U{1) symmetry of the background. 
Hence, if a consistent theory of quantum gravity can be defined in NHEK with these boundary 
conditions then it must be a chiral CFT. There has been considerable interest in extending the 
Kerr-CFT conjecture, and entropy calculation, to other extremal black holes [3]. 

The GHSS boundary conditions are unusual in two respects. First, they specify the rate at 
which components of hfj,u (the deviation of the metric from the NHEK geometry) should behave 
asymptotically. We shall refer to these as the "fall-ofF' conditions. Most components decay rel- 
ative to the background but some are allowed to be 0(1) relative to the background. Secondly, 
GHSS impose a supplementary boundary condition, namely that the energy (the conserved charge 
associated with the generator Lq of SL(2,R)) should vanish. 

One motivation for this paper is that the GHSS fall-off conditions are motivated entirely by 
considerations of the asymptotic symmetry group. However, boundary conditions are also required 
for classical physics to be predictable from initial data in a non-globally hyperbolic spacetime such 
as NHEK (or anti-de Sitter). It is not clear whether these boundary conditions will be compatible 
with the unusual GHSS boundary conditions. Indeed, it is not even clear whether the GHSS 
boundary conditions allow propagating gravitational degrees of freedom, or whether they lead to 
physics similar to Einstein gravity in AdSs, where non-trivial physics is associated with large gauge 
transformations (i.e., non-trivial elements of the asymptotic symmetry group) and black holes that 
are locally, but not globally, gauge [4 . We shall investigate these issues by studying linearized 
gravitational perturbations of NHEK. 

Another motivation for studying perturbations of NHEK is associated with positivity of the 
energy. The GHSS "zero energy" condition arises from the desire to consider only the ground 
states corresponding to an extreme Kerr black hole, rather non-extremal excitations. However, this 
presupposes that the energy must be non-negative. The NHEK geometry possesses an ergoregion, 
inherited from the ergoregion of the Kerr black hole. It is well-known that, in the presence of an 
ergoregion, one can construct initial data for test matter fields for which the energy of these fields 
is negative [5^ For a Kerr black hole, this is not a problem because the positive energy theorem 
[6] ensures that the total energy of the spacetime (black hole plus matter) is non-negative. This 
is a non-trivial result, which may not extend to NHEKJ^ Furthermore, in a spacetime with an 
ergoregion but no event horizon, e.g. NHEK (adopting the global perspective), if one imposes 
boundary conditions such that there is no energy in the matter fields entering from infinity, then 
the total energy of these fields can only decrease. If it is initially negative then it will become more 
negative, suggesting an instability [5] . 

It should be noted that the issue of NHEK stability is subtle: BH pointed out that the singularity 
theorems imply that there exist small perturbations of NHEK that will lead to the formation of 
a singularity. In this sense, NHEK is unstable. However, as BH also observed, such a singularity 

^ If one wanted to prove such a theorem using spinorial methods then NHEK would have to admit a spinor field 
covariantly constant with respect to some connection. As far as we know, no such spinor field has been constructed. 
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might be hidden inside a tiny black holej^ If this has positive mass then there would not be a 
problem. However, if the energy is negative, or the singularity is naked, then it would be difficult 
to make sense of NHEK. 

The NHEK geometry shares many similarities with AdS^: indeed, it is foliated by warped 
^^5*3 submanifolds, which have been discussed extensively in recent work on topologically massive 
gravity (TMG) [7j. In TMG, there are propagating gravitational degrees of freedom but some of 
these turn out to have negative energy, signaling a potential instability of AdS^ [8]. In the chiral 
limit, the propagating modes are eliminated by boundary conditions at infinity [H [9], leaving only 
pure gauge modes and BTZ black holes, just as in Einstein gravity . Away from the chiral limit, 
AdS3 is unstable but there exist warped AdSs solutions that might provide an alternative ground 
state [To]. The stability of some of these has been investigated recently [Tl]. Again, there are 
propagating modes with negative energy but these are excluded by boundary conditions. 

We now describe the approach we shall take. NHEK is a type D vacuum spacetime so one 
can obtain decoupled equations describing gravitational perturbations using Teukolsky's method 
|12| [T3] . The Teukolsky equation turns out to be very similar to the equation governing a massless 
scalar field in NHEK, which was discussed by BH, and the qualitative features of our solutions 
closely resemble theirs. 

By expanding in (spin-weighted, spheroidal) harmonics on the of the NHEK geometry, 
we reduce the Teukolsky equation to the equation of a charged massive scalar in AdS2 with a 
homogeneous electric field. This equation can be solved in terms of hypergeometric functions. 
Depending on the labels {l,m) of the spheroidal harmonics, the solutions either grow or decay as 
powers of the AdS2 radial coordinate, or they are oscillating at infinity. In the former case, the 
natural "normalizable" boundary conditions lead to quantized frequencies: we shall refer to these 
as normal modes. These modes fill out highest-weight representations of a Virasoro algebra which 
extends the SL{2,R) isometry group of AdS2, indeed such modes have been obtained previously 
in the context of a charged scalar in AdS2 with electric field A particularly important set of 
normalizable modes are those arising from axisymmetric (m = 0) perturbations of NHEK. 

The other set of modes are those that oscillate at infinity. Following BH, we refer to these 
as traveling waves. These modes typically have large m for given I: \m\ ~ I. From the AdS2 
perspective, these correspond to modes that have complex weight with respect to the generator 
Lo of SL(2, R) and so would not normally be considered. However, in NHEK it would be very 
restrictive to discard these modes since that would correspond to a restriction on the allowed values 
of (/, m). Even if such a restriction were imposed at the linearized level, it would be violated at the 
nonlinear level through interactions between modes. 

The traveling waves carry energy and angular momentum to infinity. BH showed that such 
modes are associated with superradiant scattering in the NHEK geometry. However, rather than 
considering scattering, we are interested in the question of what happens to localized initial data. 
We therefore impose purely outgoing boundary conditions at infinity. We find that the modes 
corresponding to traveling waves become exponentially damped, i.e., they are quasinormal modes 

■^The same might be true in AdSd for d > 4: AdSd is like a confining box, and a small gravitational perturbation 
in a box might be expected to evolve ergodically. If so, eventually sufficient energy will be concentrated into a small 
enough region to produce a tiny black hole. We thank G. Horowitz for discussion of this point. 
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of NHEK, describing the decay of a small perturbation via radiation to infinity. Therefore NHEK 
is stable against linearized gravitational perturbations. The reason that the above argument for 
instability based on the energy in matter (or linearized gravitational) fields fails is that some 
outgoing waves carry negative energy to infinity. Hence the energy flux through infinity need not 
be positive and so the energy need not decrease with time. 

So far, our discussion of gravitational perturbations has been based entirely on the Teukolsky 
equation. However, in order to calculate the energy, or discuss fall-off conditions on the metric, 
we need to know the perturbed metric tensor rather than just the Teukolsky scalars. Fortunately, 
there exists a method for determining the metric perturbation in terms of a scalar potential, called 
the Hertz potential |15j-jl9j. This satisfies an equation closely related to the Teukolsky equation. 
Using this, we obtain explicit results for the form of the metric perturbation. 

We find that most (but not quite all) normal modes satisfy the GHSS fall-off conditions but trav- 
eling waves violate these conditions. Although one can construct localized wavepackets involving 
the latter, they will eventually propagate to infinity and violate the fall-off conditions. Therefore, 
at the linearized level, they should be excluded, leaving just the normal modes. 

Next, we consider the energy of the normal modes. To warm-up, we start by considering a 
massless scalar field. We are able to show that an arbitrary superposition of normal modes has 
positive energy. Then we turn to gravitational perturbations. We define the energy of the latter 
in the usual way using the Landau-Lifshitz "pseudotensor" . Since the metric perturbation involves 
second derivatives of the Hertz potential, the energy involves an integral of a complicated quantity 
sixth order in derivatives. Nevertheless, using a combination of analytical and numerical methods, 
we find that the energy of gravitational normal modes is positive, thus supporting the validity of 
the GHSS zero-energy condition. 

This positive energy result is satisfying but the exclusion of the traveling waves is worrying. 
First, it is worrying that we can construct initial data that satisfy the fall-off conditions, but violate 
these conditions when evolved. It suggests that the initial value problem, at least for linearized 
fields, may not be well-posed. Furthermore, if one goes beyond linearized theory then interactions 
between modes will excite traveling waves even if they are not present initially]^ So one might worry 
about well-posedness of the nonlinear theory too. It is possible that these problems are cured by 
backreaction, i.e, going beyond the linearized approximation. We shall discuss this further at the 
end of the paper. 

This paper is organized as follows. In section [2| we derive and solve the Teukolsky equation in 
the NHEK background, obtaining the spectrum of normal, and quasinormal modes. In section[3]we 
introduce the Hertz potential and use it to obtain the explicit form of linearized perturbations. We 
compare the asymptotic behaviour of these with the GHSS boundary conditions. We then calculate 
the energy of scalar field and gravitational normal modes. Finally, section |4] discusses how going 
beyond the linearized approximation may solve some of the problems just discussed. 

Note added. As this work was nearing completion, we learned that another group is exploring 
similar issues |2U] . 

^ The only way to escape this conclusion is to consider only axisymmetric (m = 0) modes, which form a consistent 
truncation of the full set of modes. 
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2 Massless fields of arbitrary spin in NHEK 



2.1 NHEK and its Newman-Penrose tetrad 

In global coordinates the NHEK metric is [1 , (we use the notation of Ref. [2] and, because we shall 
employ the Newman-Penrose formalism, a negative signature metric) 

ds^ = -2Gm'^(e) -(1 + r2)dt2 + + de^ + K^ieW + rdtf , (2.1) 

\ 1 -|- J 

with 

n\e)^\{i + cos'e), a(0) = -4^, gj = g^mI,,^^m\ (2.2) 

2 1 -|- cos^ 

Surfaces of constant 9 are warped AdS^ geometries, i.e., a circle fibred over ^^5*2 with warping 
parameter Is?{6). The isometry group is SL{2, R) x U{1). BH showed that the solution is geodesi- 
cally complete, with timelike infinities at r = ±00. There is an ergoregion (where d/dt is spacelike) 
which extends to r = itcxD. 

In the next subsection we study perturbations in the NHEK using the Teukolsky formulation. 
For that we need the Newman-Penrose (NP) tetrad, spin coefficients and directional derivatives. In 
Appendix |A] we obtain the shear-free null geodesies of this background and use them to construct 
the associated NP null tetrad jJT], e^^) = e(2) = n, e(3) = m, 6(4) = m*, where (coordinates are 
listed in the order {t, r, 6, cj)}) 

m^ = ^ (0,0,l,iA-\e)) , (2.3) 

and e^"^) = r/^'^^^^'^ef^jj^ with non-vanishing symmetric rj^"-^^''^ = 'n{a){b) given by rj^^^^'^^ = — 7y(3)(4) _ -|^_ 
This NP tetrad satisfies the normalization and orthogonality conditions (A. 16), and the null vector 
£ is tangent to affinely parametrized geodesies: i^Vf^iu = 0. 

The unperturbed Weyl scalars in the NHEK geometry are computed using (A.20), yielding 

\I/o = = = ^1^4 = , 

yii2 = -[M^{l-icos9f]~\ (2.4) 
The first line confirms that this solution is indeed Petrov type D. 



2.2 Teukolsky master equation 

Teukolsky has shown how, for type D spacetimes, one can use the NP formalism to derive a 
system of decoupled equations, that furthermore separate into an angular and radial part, for the 
perturbations of several NP scalars [121 [13]. For gravitational perturbations, the relevant quantities 
are the perturbed Weyl scalars ^'q^^ (spin s = +2) and ^'4"'^^ {s = —2); the complex NP scalars (/>o,i 
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for spin s = ±1 Maxwell perturbations; the Weyl fermionic scalars xo,i for massless spin s = zt^ 
perturbations; and the scalar field ^ for massless spin s = perturbations. Teukolsky's master 
equation encompasses all of these cases [13] . 

Using the NP quantities listed in appendix |A] we find that the Teukolsky master equation for 
spin s field perturbations 'i/^^^ in the NHEK geometry is 



1 



2r 



(l+r2) * (l+r2, 



) ^ \ 1 + 4 sm"' 6 j ^ 



sm ( 



doismOdo^^'A -2s 



fl+r2 



-2^ (jV^ + + 4 cos ^1 ^-/-^^'^ + cot' ^ - ^^'^ = • (2-5) 

\(l+r^j sm 6* 2 J ^' 

We have allowed for the possibility of a source term on the RHS (see Appendix [C]) . The relation be- 
tween the nomenclature used here and the orig inal notation of Teukolsky [13] is {^'(^), ^'(i), ^'(1/2)} = 
{^'o, </>0; Xo} and {T(2), T(i/2)} = {^o, Jq^T^q} fo^ positive spin. For negative spin the map is 
{\I/(-2)^ vj/{-i/2)| ^ {(-^2)^ ^4, (-^2)^ <A2, (-^2)^ Xi}- Here, the powers of the unper- 

turbed Weyl scalar ^'2 are those that allow for the separation of the master equation, when we 
further assume an ansatz for the perturbation that is a radial function times the spin-weighted 



spheroidal harmonic; see (2.6). For the source term one has the map {T(_2)) ^(-^1)1 ^(-1/2)} 
{T4, J2,T^^}. These relations are summarized in Table [T| 
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(-^2)^ 02 


00 
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Xo 
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-2 
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-1 
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1 
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1 
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2(-'J'2)3T4 


2To 


(-^2)3 J2 


Jo 


(-*2)3T^l 


T 





Table 1: Teukolsky fields spin s and source terms for the master equation (]2.5|) 



2.3 Separation of variables 

We shall solve the Teukolsky equation in the NHEK geometry by separation of variables. Assuming 



""*e™^^;i(r)5;;V) (-^2)-^ , s<0, 



(2.6) 



equation (2.5) separates into an angular and radial equations. The angular equation is 

\2" 



sin^ de 



1 d f . „d 
sm 



(C cosey-2sC COS9 + S + A 



(s) [m + s cos ( 



Im 



sm 



(2.7) 
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for C = m/2 and where is the separation constant. Its eigenfunctions are the spin- weighted 
spheroidal harmonics e,'^^'^ s[^{9) (the nomenclature usually includes an appropriate normalization 
factor; see e.g., [22]), with positive integer / specifying the number of zeros, i — max{|m|, of 

is) 

the eigenf unction. The associated eigenvalues can be computed numerically with very good 
accuracy and are specified by s,l,m subject to the regularity constraints that —l<m<l must be 
an integer and / > The transformation 6 ^ ir — 6 can be used to show that 



A 



Im 



A 



A 



Im 



We also note that, to leading order in C, A, 
\m\ <C I. 



m 



{I - s){l + s + 1) + 0{C). This is useful when 



Equation (2.7) represents the most standard way to write the spin- weighted spheroidal harmonic 



(s) 

equation. However, it will be convenient here to work with shifted eigenvalues A, defined by 



A 



Im 



+ + s 



The advantage of using these quantities is that they have the symmetry 



A 



i-s) 
Im 



A 



is) 
Im ' 



(2.9) 



(2.10) 



Notice that in the Kerr background with mass M and angular velocity the angular equation 

= auj, where a = 2Mr^Q}{ is Kerr's rotation 



for spin s perturbations is also (2.7) but with Ck 



parameter and oj the wave's frequency in this geometry. As observed in [1 , in the near-horizon limit 
of extreme Kerr, all finite frequencies uj in the NHEK throat correspond to the single frequency 
to = mO|^* = ^ in the extreme Kerr geometry. This uj corresponds precisely to the marginally 
unstable superradiant frequency, and in the NH limit one finds CKcrr = Muj C = m/2. 
Writing for any spin. 



(2.11) 



we find that the radial equation associated with ( |2.5| ) can be written also in a unified way as 

{oj + qr)'^ 



dr 



^["'^ frl = 



(2.12) 



with 



q = m — IS , 

,' = <f + ^2=^\2+s-2^sm-^^ 



(2.13) 



This is exactly the equation for a charged massive scalar field in AdS2 with a homogeneous electric 
field: take the AdS2 metric in global coordinates, 

dr2 



dsi = {l+ r'^)dt^ 
and the electric field to arise from the potential 

A = rdt. 



l+r' 



;dr\ 



(2.14) 



(2.15) 
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Figure 1: Values of rj , defined in (2.20), for \s\ =2, and a) / = 2 and b) / = 3. 
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Figure 2: Values of r/^, defined in (2.20), for |s| = 2, and a) / = 4 and b) I = 16. 



Define the covariant derivative for a field of charge q as 

P = V - iqA, 



(2.16) 



where V is the Levi-Civita connection in AdS2- The equation for a charged scalar field ^{t, r) with 
mass /i is then 

(2.17) 



V'^<^ + = 0. 



Assuming 



$(t,r) 



-ILOt 



$(r). 



(2.18) 



the equation of motion reduces to (2.12). Therefore, a general spin s perturbation with angular 



momentum m in NHEK obeys the wave equation for a massive charged scalar field in AdS2 with a 
homogeneous electric field. However, note that the charge q is complex, as is the squared mass 
although /i^ — q^ is real. The problem of a massive charge scalar field in AdS2 with homogeneous 
electric field was studied in Ref. [13!, where solutions corresponding to highest weight representa- 
tions of a Virasoro algebra extending SL{2, R) were obtained. We shall recover the same solutions 
in the next section. 
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2.4 Solving the radial equation 



Asymptotically, the solutions of (2.121 behave as 

$(r) ~ |rr^/2±W^ (2.19) 

where 

7] = Vl+4(/i2-g2) = ^1 + 4AJ^, Im(77) > 0. (2.20) 

Note that r](s,l,m) = r](—s,l,m) = r]{s,l, —m). We can now see that the modes can exhibit 
qualitatively different behaviour, depending on the value of (Z,m), as first noticed by BH (for 
s = 0). Some modes have real rj and others have imaginary rj. For example, axisymmetric modes 
(m = 0), have, for general s, 

?? = 2/ + l (m = 0). (2.21) 

i.e., such modes exhibit power-law behaviour at infinity. However, for certain other modes, specifi- 
cally those with |m| ~ /, rj is imaginary and hence the solutions oscillate at infinity. In Figs. [Tjand 
[2] we show how tj'^ depends on m for gravitational perturbations with some different values of I. 

It is interesting to ask which modes have the smallest real value for r] since these will give the 
normal modes that decay most slowly at infinity. For gravitational perturbations (|s| = 2) we have 
calculated rj for all {Um) with / < 30 and find that the mode with the smallest real value for r/ 
occurs for I = 4, \m\ = 3, which gives r] = 2.74. 



Equation (2.12| can be solved exactly. This is not a surprise since in the Kerr geometry, 
Teukolsky and Press [23] found that the corresponding Teukolsky radial equation can also be 
analytically solved in the particular case where we have extreme Kerr and a wave frequency that 



saturates the superradiant bound, ui = m^l^*. As discussed after (2.10), all frequencies in NHEK 
correspond to the single superradiant threshold frequency in the extreme Kerr. So we indeed expect 
this property for the radial equation in NHEK. 



Introducing the new radial coordinate. 



and redefining the radial wavefunction as 



1 - ir) , (2.22) 



^Sl(0 = ^"(l-^)^i^, with a^^{u;-iq), i3^^{u; + iq), (2.23) 



the radial equation (2.12) can be rewritten as 



z{l - z)d'lF + [2a - 2 (a + (5) z] d^F - [(a + /3 + l)(a + /3 - 2) + (g^ - /i^ + 2)] F = . (2.24) 

This wave equation is a standard hypergeometric equation |24] . z(l — z)d'^F+[c— {a + b + l)z]dzF — 
abF = 0, with 

a= ^(l + r/ + 2w) , 6= ^(1 -r? + 2w) , c = 1 + w - , (2.25) 
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and hence the most general solution in the neighborhood of z = is 

= - zfF{a, b, c, z) + Bz'^+^-^{l - zfF{a - c + l,b - c + 1,2 - c, z). (2.26) 

We render this function single valued in the complex z plane by taking branch cuts to run from 
— oo to and from 1 to +oo, corresponding to taking | arg(2;)| < tt, | arg(l — z)\ < vr. Note that the 
branch cuts do not intersect the line Re(z) = 1/2, which corresponds to real r. 

2.5 Boundary conditions 

The above solution of the radial equation is regular for all finite r. Using standard properties of 
the hypergeometric function, we find that it exhibits the following behaviour as r — > ±00: 



2 



(2.27) 



where 



= A-j:^ - Be^- r(2-c) 



T{b)T{c-a) r(6-c+l)r(l-a)' 
n± _ A r(c) p^ifcCTc r(2 - c) 

^ ~^T{a)T{c-b) r(a-c + l)r(l-6)- ^^■^^> 

The boundary conditions now depend on whether rj is real or imaginary. 

2.5.1 Normal modes 

Assume that rj is real. In this case, we impose normalizable boundary conditions, corresponding 
to demanding that = = 0, a pair of simultaneous equations for A, B. Non-zero solutions 
exist only if the determinant of this system vanishes. Using r(2;)r(l — z) = tt/ sm^nz), this gives 

- = (2 29) 

r(a)r(i-6)r(c-6)r(a-c + i) ' ^' ' 

This imposes a quantization condition on oj, corresponding to the two solutions a = —n and 
1 — b = —n where n = 0, 1, 2, . . .j^ The former solution gives 

a; = -(n+l/2 + r//2), n = 0,1,2,..., 5 = 0, (2.30) 

and the latter gives 

w = n + 1/2 + r//2, n = 0,1,2,..., ^ = 0. (2.31) 
We can summarize the normal mode spectrum as 

(J = ±(n+ 1/2 + r//2) , n = 0,l,2,.... (2.32) 



■'At first sight, condition (2.291 could also be satisfied if we imposed c = 1, i.e., u) = iq. However, a more careful 



analysis rules out this possibility because for c = 1, (2.261 is not a solution of the problem: one must allow for a 



logarithmic dependence in the second part. Redoing the analysis with the appropriate regular radial solution for this 
special case [21], we conclude that nothing physically special occurs for c = 1. 
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This is precisely the spectrum of normal modes found for a massive charged scalar in AdS2 with a 
homogeneous electric field in Ref. jl4j . 

Note that we have allowed uj to be positive or negative. This is because the Teukolsky equation 
for s 7^ is not invariant under complex conjugation, so negative frequency solutions are not 
simply related to positive frequency solutions by complex conjugation, they have to be considered 
separately. The two possible signs correspond to the two different helicities of the field. The radial 
equation is invariant under uj —>■ —u>, r ^ —r hence ^[^(.^jC'') ^ ^imuji~''^)- 

For n = 0, the positive frequency solution of the radial equation is 

$[^(,=0) ^ _ ^)-(l+r,)/4-V2^ (2.33) 

and the negative frequency solution is obtained by r — > — r, i.e., z ^ 1 — z. The solutions with 
positive n are related to these n = solutions by multiplication by a polynomial of degree n in z. 



2.5.2 Traveling w^aves 

Now consider the case of imaginary r]. Define fj > hy r] = ifj. The radial function oscillates at 
infinity, corresponding to incoming or outgoing waves (see Appendix [B] for details). Rather than 
considering scattering in NHEK, we shall impose boundary conditions corresponding to purely 
outgoing waves at infinity, which will discretize the frequency uj and render it complex. A solution 
with positive imaginary part corresponds to an instability, and a solution with negative imaginary 
part is a quasinormal mode. 

As discussed by BH, there are two inequivalent notions of "outgoing" that one can use in NHEK 
because the phase velocity and group velocity of wavepackets need not have the same sign, e.g. for 
positive UJ and m, the group and phase velocities have the same sign at r = +oo but opposite sign 
at r = — oo (see Table [3]). Physical boundary conditions correspond to the notion of "outgoing" 
defined using the group velocity. However, it is easier to analyze the case of outgoing phase, so we 
shall consider this case first. 

Assume that Re(a;) > 0. Then the solutions with outgoing phase at r ^ ±00 are the solutions 
with = 0. This leads to the quantization condition 

(1 — c)tt 

^ ' -0, (2.34) 



r(6)r(l - a)T{c - a)T{b - c + 1) 



with solution 1 — a = —n, n = 0, 1, 2, . . . (6 = —n is inconsistent with Re(w) > 0), which gives 
u! = n + 1/2 — ifi/2. Repeating the exercise for Re(u;) < requires = 0, and leads to u; = 
— (n + 1/2) — ifj/l. We can summarize the result as 

UJ = n + l/2 -ifi/2, n£Z (2.35) 

The imaginary part is negative, hence these are quasinormal modes. This is a little surprising. BH 
pointed out that the energy flux (for positive frequency modes) has the same sign as the phase 
velocity. Hence outgoing phase should correspond to outgoing energy at infinity. As discussed in 
the introduction, this is precisely the situation in which one expects an instability associated with 
the negative energy in matter fields within the ergoregion becoming increasingly negative. We have 
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found that outgoing phase leads to stable quasinormal modes rather than an instability. However, 
these boundary conditions are unphysical: we are arranging that an initial wavepacket (at finite r) 
composed of modes with positive a;, m does not propagate to r = — oo by sending in an appropriate 
(finely tuned) wavepacket from r = — oo to scatter with it in such a way as to produce only a 
wavepacket propagating to r = +00. This is analogous to boundary conditions for a Kerr black 
hole in which one arranges that initial data leads to no waves crossing the future horizon by sending 
in appropriate waves from the past horizon. Presumably, the fine-tuning is the reason that we do 
not see an instability here. 

Now consider the physical boundary conditions corresponding to "outgoing" defined with re- 
spect to the group velocity. Assume that Re(u;) > and m > 0. BH showed that, under these 
conditions, the phase and group velocities have the same sign for r ^ 00 but opposite sign for 
r — > —00. Hence the boundary conditions that we need are C"*" = D~ = 0. In fact, the same holds 
for Re(a;) < and m > (see Appendix [B|) . Using the identity r(z)r(l — z) = it/ sin(7rz), we find 
that the quantization condition is 



sin(7r6) sin[7r(c — a)]e 



sin(7ra) sin[7r(c — b)]e'''"'^, 



which give^ 



1 ^ 1 

LJ = n-\ log 

2 2iT ^ 



cosh[7r(?7/2 -|- m)] 
cosh[7r(r)/2 — m)] 



n e Z, 



(m > 0) 



(2.36) 



(2.37) 



where we have specialized to scalar field (s = 0) or gravitational (±2) perturbations for simplicity. 
Repeating the analysis for m < requires = = 0. The general result is 



cosh[7r(7//2 -|- \m\ 
cosh[7r(7//2 — \m\ 



n G Z 



(2.38) 



We see that lm{uj) < hence these are stable quasinormal modes. So NHEK is stable against 
linearized gravitational (and scalar field) perturbations. 



3 Metric perturbations 
3.1 Hertz potentials 

The NP scalar perturbations '^/^'^'^ are useful because they are invariant under infinitesimal diffeo- 
morphisms and under rotations of the NP tetrad. Many physically interesting quantities can be 
computed directly from the knowledge of these NP fields |12| [T3\ [2T] . However, in some problems as 
is our case, we really need to know the perturbations of the metric itself, h^^, or the perturbations 
of the Maxwell or Weyl fermionic vector fields, respectively and x^l■ Cohen and Kegeles |15|I17|. 
and Chrzanowski [16] have proposed a unique map that provides the h^,^, or Xfi perturbations 
given the so-called Hertz potential ^'[^^ (see a good discussion also in [19j). Wald proved Cohen- 
Kegeles— Chrzanowski 's results [18]. See Appendix [C] for a detailed discussion of these works. The 

solution corresponding to c taking integer values is ruled out for the reason discussed in footnote 4 
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main conclusion is that the Hertz potential also obeys a pair of decoupled equations, again one for 



positive and the other for negative s. These are written in equations (C.18) and (C.19). 

For gravitational perturbations, this method yields the metric perturbation in a particular 
gauge: the ingoing (outgoing) radiation gauge IRG (ORG), specified by the conditions 

i^h^u = g^'Ku = (IRG) , V = a^'Ky = (ORG) . (3.39) 

At first sight, these gauge conditions appear overdetermined but it has been shown that, for pertur- 
bations of a type II vacuum spacetime, there is a residual gauge freedom that allows one to impose 
the IRG provided that l^i^T^y = where £ is the repeated principal null direction and T^^ the 
stress-tensor of any matter perturbation present j25j . Similarly, for type D one can impose either 
the IRG or the ORG (if n^^-nyT^y = 0). The spin of the Hertz potential corresponds to these two 
different gauges: the metric perturbation in the IRG (ORG) is obtained from the Hertz potential 
with s = —2 (s = +2). The two Hertz potentials contain exactly the same physical information, so 
one need only work with one of them. 

For vacuum type D spacetimes, the Hertz potential itself satisfies a master equation. For the 
Kerr solution, this master equation turns out to be exactly the same as for the original NP scalars 
"^^^\ equation (2.5 1, with no source term on the RHS [261. We have checked that the same is true 



for NHEK. More concretely, ^''^ = {^'^ ^\ ^''^''i the Hertz potentials conjugate to 

the positive spin Teukolsky perturbations \j/(i/2)| ]-,^^ satisfy exactly the same master 

equation (2.5) as (—^2)"^^^^^ for negative spin. Similarly, ^jf^ = ^jj^^^} are the 

Hertz potentials conjugate to the negative spin Teukolsky perturbations vj/(-i/2)| 
but the positive spin Hertz potential (— ^'2) ^ obeys the same master equation as for 
positive spin. In short, the Hertz potential obeys the same master equation as its conjugated 
Teukolsky field but with spin sign traded. This relation is better clarified if we use Tables [1] and |2] 









(-^2)-^< 


Hn 


(-M/2)-^.I;W 


Hn 


(-^2)-HfJ/^) 




s 




-2 


2 


-1 


1 


1 
2 


1 
2 






is) 

Table 2: Spin s Hertz fields ^'jj that satisfy the master equation (2.5) with no source term. 



Assuming perturbations for the Hertz potentials of the form 



(s) 
H 



-iuji^irri' 



(3.40) 



where Rimuj{T) further satisfies (2.11 ), equation (|2.5[) separates into an angular and radial equations. 



The angular equation is (2.7) with C = m/2, and the radial equation is (2.12). Its solution is given 



by (2.26). 



As stated above, given the Hertz potential for the gravitational field there is a unique map 
between it and the metric perturbations jl5| [TH| [T71 [TS] . A similar map exists between the spin 
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s = ibl,ibl/2 Hertz potentials and the Maxwell and Weyl fermionic vector perturbations, but 
we leave the discussion of these cases to Appendix [Cj In the ingoing radiation gauge the metric 
perturbation in NP notation is given by (see Appendix O 



h 



IRG 



^(^m^) [(D + 3e + e - p + + 4/3 + 3r) + ((5 + 3/3 - a - r - lt){D + 4e + 3p)] 



+ 3/3 + a - r)((5 + 4/? + 3r) - m^m^(Z) + 3e - e -/))(£> + 4e + 3p) |*h + c.c. , 

(3.41) 

and a similar correspondence exists between the Hertz potential and the metric perturbations hPj}''^ 



in the outgoing radiation gauge. (See the second relation of (C.IO).) One can check that (3.41) 



indeed satisfies the linearized Einstein's equations for a traceless metric perturbation: 



(3.42) 



3.2 Behaviour of solutions 



The basis vector fields I and n are globally well-defined. However, the vector field m is singular at 
9 = 0,TT. Nevertheless, one can check that angular dependence of the Hertz potential contains a 
sufficiently high power of sin^ to ensure that the above metric perturbation is smooth at 6 = 0,it. 

(±2) A 1 A 1 

The asymptotic behaviour of the Hertz potential can be obtained using (2.11) and (2.19). 

Use of (3.41) yields then for the asymptotic /i^jy behaviour (rows and columns follow the order: 

/ 0(1 



{t,r,e,<p}) 



Oil: 



\ 



Oil 
Oil) 



(3.43) 



where rj is given by (2.20). Exactly the same result is obtained in the outgoing radiation gauge. In 
(3.43) we have not imposed any boundary condition. These were discussed in subsection 2.5 e.g., 



for rf > 0, the lower sign would correspond to normal modes. 

We shall now compare the above asymptotic behaviour of metric perturbations with the GHSS 
fall-off conditions. The tr and t9 components are the most restrictive. For these to satisfy the fall- 
off conditions, r] must be real, so traveling waves are excluded, we must use normalizable boundary 
conditions (i.e. the lower sign choice) and we need r] >3. Recall that there are normal modes with 
7] = 2.74, so it appears that the GHSS fall-off conditions exclude some of the normal modesj^ 

As emphasized in the introduction, at the nonlinear level, we expect that interactions will lead 
to modes corresponding to traveling waves (r/^ < 0) being excited, which would lead to a violation of 
the GHSS fall-off conditions. The only modes that escape this conclusion are the axisymmetric ones 
(which have with r] = 21 + 1), which always obey the GHSS boundary conditions. Axisymmetric 
modes form a consistent truncation of the full set of modes in the sense that linearized axisymmetric 
modes will not excite non-axisymmetric modes at next order in perturbation theory. 



It is conceivable that a gauge transformation could be used to bring a mode violating the fall-off conditions to 
one that satisfies these conditions but this seems unlikely, especially for traveling waves. 
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3.3 The energy 

3.3.1 Massless scalar field 

We want to compute the energy associated with the gravitational perturbations that we found in 
the previous subsection. Since this wih involve a rather lengthy calculation, we shall start with the 
conceptually simpler case of a massless complex scalar field: 

□ ^> = 0. (3.44) 

The canonical energy momentum tensor is given by 

T^u = V(^W,)r - ^^^.V«$V"$*. (3.45) 

Let S be a spacelike hypersurface with future-directed unit normal n^. Then, given any Killing 
vector we can define the associated conserved charge 

= [ (fx^T'^'^n^^,, (3.46) 

where = g^y — n^riy is the induced metric on S. We shall choose S to be a surface of constant 
t in the NHEK geometry. The conserved charges of interest are the energy for ^ = d/dt, and 
the angular momentum J , for ^ = —d/d(j). (The latter is the U{1) charge of GHSS.) Written out 
explicitly, these are 



£ = M'^ d(l) de dr 

Jo Jo J~oo 

J = M'^ I d(j) d9 dr 

Jo Jo J~oo 



l+r2 ' ^ ' yi-'^ I ' 1--^ I ' \H0)^ l+r2^ 
^ ^2r \d^<^\^ - {dt^d^^* + dt^*d^^)] . (3.47) 



1 + r 



In the energy integrand, all the terms except the last are manifestly positive. The last is proportional 
to = gtt and thus is positive only outside the ergosphere where 1 + — r^A(0)2 > 0. The 
energy can thus be negative (for a rigorous proof of this, see Ref. [5j). 

Consider a general superposition of normalizable modes (recall that n is defined by the frequency 



quantization (2.32)): 



^(^) = Yl i'^nlm'^nlmix) + bnlm'^ nlmix)* ) ■ (3.48) 
nlni 

First, we shall show that the conserved charges associated with such a solution can be decomposed 
into a sum of conserved charges of the individual modes. 

A charge integral can be regarded as defining a (typically indefinite) norm on the space 

of solutions of the wave equation. Given a norm ||, it can be "polarized" to obtain a Hermitian 
scalar product (,): the real and imaginary parts of {u,v) are given by {\u + v\ — \u\ — \v\)/2 and 
{\iu + v\ — \u\ — \v\)/2 respectively. In our case, polarizing the charge integral defines a scalar 
product (<I>i, <I>2)g) antilinear in $1 and linear in $2- Since the norm is conserved, so will be the 
scalar product. Note that (<I>,<I>)5 = 
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We shall now argue that modes with different {nlm) are orthogonal with respect to this scalar 
product. The scalar product has the form 

($1, $2)5 = ^ d'^xQ^"'{x)^^<^l^u'^>2, (3.49) 

where Q^^ is preserved by any Killing vector field that commutes with ^. Now let r] be such a 
Killing field. We can then write 

($1, -i£^$2)c - {-iCrj^uhh = Cr,{Q^''{x)d^^\d^^2). (3.50) 

If the RHS vanishes then this shows that —iCrj is self-adjoint with respect to the this scalar product. 
For NHEK, we take ^ = d/dt or d/d(j). Taking r/ = d/dt^ the RHS vanishes because the scalar 
product is conserved, and hence independent of t. Taking ry = d/d(j), the RHS vanishes because it 
is a total derivative on S. It follows that modes with different uj or diffferent m will be orthogonal 



with respect to this scalar product. Hence, in calculating the charge associated with (3.48), there 



are no cross-terms in the charge arising from modes with different w or m (in particular, there are 



no cross-terms between the positive and negative frequency parts of (3.48)). 

Now consider the /-dependence. Since I is not associated with a Killing symmetry of the 
background, we cannot use the above argument. Instead, for separable solutions, the angular 



dependence will be given by (2.7) with s = 0. This equation is self-adjoint, so two solutions with 



(s) 

different values of will be orthogonal with respect to the measure s'mO, i.e.. 



rdesmesl'^l{e)sl^l{er oc 6i,i,. (3.51) 

Jo 



Fortunately, it turns out that sin^ is precisely the measure that arises in the scalar products 
associated with the energy and angular momentum. 

From these results, we see that no cross-terms between modes with different (nlm) contribute 



to the energy and angular momentum. Substituting (3.48) into (3.46) for S, = dt gives the energy 
as a sum over contributions from individual modes: 

{\anlm\ + \bnlm\ ), (3.52) 

nlm 

where 



4vrM2u;„,„/ d0sin^|5S(0)|W dr|iiil(r)p5^^=?^ . (3.53) 



Jo J-oo i + r 

Note that Snim is manifestly positive only when m = 0. However, we have evaluated the radial 
integral above for many cases, namely for < / < 10, — / < m < / and < n < 10. In all these 
cases, it is positive. Hence, for a massless complex scalar field in the NHEK geometry, the energy 
of an arbitrary superposition of normalizable modes is positive. 
The angular momentum can be similarly decomposed: 

J'nlm{\^nlm\ l^nZml ); (3.54) 

nlm 

where we find the simple result 

J^nlm 



^nlm ^nlm 



(3.55) 
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3.3.2 Gravitational perturbations 



The energy of gravitational perturbations is calculated from the Landau-Lifshitz "pseudotensor" 
defined as follows. Consider metric perturbations h^i, around NHEK up to second order in the 
amplitude, 

9^^u = 9^.u + V = 9^^r. + h'jjj + /ig) +0{h^), (3.56) 
The linearized Einstein equation i^ 

G«[/iW]=0. (3.57) 
At second order, the Einstein equation relates terms linear in h^"^^ to terms quadratic in h^^^: 



(3.58) 



where the RHS is quadratic in h^^\ Written out explicitly, for traceless perturbations it reads 
(here, we use the notation 



"fMl/ ll'lJ.U ) 



8-kGT, 



¥7 



(3.59) 



We now define the conserved charges (^^[/i^^-*] associated with the first order perturbation exactly 
as in (3.46) with = d/dt, —d/dcj) giving the energy and angular momentum respectively. 

Recall that /i^^-* is related to the Hertz potential by equation (3.41), which is second order in 
derivatives. It follows that the conserved charges are given by integrals of quantities that are sixth 
order in derivatives. Hence calculating these charges involves very lengthy calculations, which we 
have performed using computer algebra. 

We consider a Hertz potential corresponding to an arbitrary superposition of normal modes: 



(3.60) 



nlm 



where s = ±2 and the superscript it refers to positive and negative frequency respectively. 

As in the scalar field, case the conserved charges can be used to define a scalar product (,)^ 
between solutions of the linearized Einstein equation. The only significant difference here is that 
the metric perturbation is real, so the scalar product will also be realj^ One can argue exactly as in 



^Written out, this takes the standard Lichnerowicz form. If we assume that /i'^^ is traceless then this equation 
reduces to (3.42 I. 

*The polarization formula is («, v) = (l/2)(|ii + u| — |m| — Ifl). We could have chosen to work with complex modes 
h^u, for which the negative frequency modes are complex conjugates of the positive frequency modes. However, then 
we would have had to take account of two different polarizations for the gravitational modes. 
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the scalar field case that Cri is anti-self-adjoint with respect to this scalar product if 77 is a Killing 
field that commutes with ^: {h,Cr^k) = {—Cr^k,h). It follows that the operator — £^ is self-adjoint 
and hence linearized metric perturbations with different o;^ or different must be orthogonal. 
Given the complexity of the charge integrals, we have not succeeded in demonstrating that 

modes with different I are orthogonal in the same way that we did for the scalar field. However, 

{'2) 

note that oo depends on Z in a very complicated way (through the eigenvalues which must 
be found numerically). Hence it seems very unlikely that modes with different / could have the 
same uj"^. Therefore the orthogonality of modes with different a;^ should ensure the orthogonality of 
modes with different I. An exception are the axisymmetric (m = 0) modes, which have uj = n + l + \ 
so modes with the same n + I have the same w. However, the axisymmetric modes are the "least 
dangerous" as far as the possibility of negative energy is concerned so we shall not worry about 
this further, and simply assume that all modes with different I will be orthogonal. 

We now turn to our calculation of the conserved charges associated with individual modes. 
These charges are most easily computed in NP tetrad, since this is the basis in which the metric 
perturbation takes the relatively simple form (3.41), although the explicit expressions for the com- 
ponents are still too long to be written here. Using Mathematica, the separated equation of motion 
can be used to reduce the integrands of the charge integrals to expressions first order in derivatives, 
which were then calculated numerically using Mathematica's NIntegrate function. 

We have calculated the energy of normal modes with / = 2, 3, 4, 5, 6 for all allowed values of m, 
n = 0,1,2,3,4,5,6 and both positive and negative frequency. In all cases it comes out positive. 
This is the main result of this section. 

The only difference between positive and negative frequency modes is the sign of J so we focus 
on the positive frequency case. The numerical value of the energy depends on the normalization 
of the Hertz potential. However, the ratio of the conserved charges, Jnim/£nim is normalization 
independent. This ratio as a function of m and n for fixed I is displayed in Figs. (|3])-(|4|. Note that 
J'ni(-m) = —Jnim and that £ni(-m) = £nim- Moreover, for any given I the ratio \ J\/£ always has 
a (non- vanishing) minimum at |m| = |s| = 2 (the ratio decreases with n but only slowly, so this 
is not apparent in the plots). The corresponding modes also exhibit special behaviour in the Kerr 
geometry: Teukolsky and Press |23] found that for a given black hole rotation and wave frequency, 
the modes whose energy is most absorbed or superradiantly amplified are precisely those with 
/ = \m\ = \s\. 



4 Discussion: second order perturbations 

We now turn to the question of what happens if we go beyond first order in perturbation theory. 
The second order metric perturbation h^"^^ is determined by solving (3.58|)r We are not going 



to attempt to solve this equation. Instead, following recent work on TMG ^27] we consider the 
conserved charges. 

So far, we have worked with conserved charges defined via bulk integrals quadratic in h^^\ 



"if we considered scalar field perturbations then the gravitational backreaction of the perturbation would be 
governed by the same equation with T^^ the scalar field stress tensor. 
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Figure 3: Ratio of the conserved charges J jE for spin |s| = 2 perturbations as a function of the 
azimuthal angular number m for a) / = 2 and b) Z = 3. We only consider normal modes, i.e., values 
of m that yield r/^ > as defined in (2.20). Data points corresponding to n = 0, 1,2 are plotted, 
with the solid line representing n = 2 and the dashed line n = 0. 



1 = 4 



1 = 6 





Figure 4: Ratio of the conserved charges J jE for spin |s| = 2 perturbations as a function of the 
azimuthal angular number m for a) / = 4 and b) ? = 6. 



However, one can also define conserved charges via boundary integrals, indeed these are the charges 
discussed by GHSS. So first we shall explain how they are related to our bulk integrals. Consider 
a 1-parameter family of exact vacuum solutions (?(A), where ^(0) = ^ is the NHEK metric. Let 
= £,'(0) and }P- = (l/2)5"(0). (/i^^) is the linearized solution arising from the linearization of 
^(A), /i*-^^ is the second order correction.) Owing to the unusual fall-off conditions, the conserved 
charge (3^(A) = (^^[^(A)] associated to a generator ^ of the asymptotic symmetry group are defined 
by integrating the following expression: 

dQi. 



Qd9'W,9W] (4.61) 



where 



/(9E 

+ ^/i""(V'^C<x + V,e^)]da;" Adx^. (4.62) 



19 



Now, our gravitational normal modes decay sufficiently fast that they give Q^[h^^\g] = 0, hence 
dQ^/dX = at A = 0. This is no surprise since we know that the energy should be quadratic in 
h^^\ Hence we have to go to next order, and calculate {l/2)d'^Q^/dX'^ at A = 0. This can be done 
by differentiating (4.611, which gives a sum of a part linear in h^'^\ equal to Q^[h^'^\g], and a part 



quadratic in h^^\ However, the normal modes decay so fast that this second part vanishes. Hence, 
to second order in A, we have that 

Q^{X) = X'Qdh^^\g]. (4.63) 
Now, assuming that ^ is a Killing field of the background, a standard manipulation \28\ [29l [30] 



based on the second-order Einstein equation (3.58) enables one to rewrite this surface integral as 



the bulk integral quadratic in h^^^ that we used in the previous section. 

One subtlety is that the NHEK geometry has two boundaries (at r = ±00). The bulk integral 
for the charge will be the sum of the two surface integrals. Hence, if the first order perturbation 
gives a non-zero conserved charge, then the second order perturbation h^"^^ must decay sufficiently 
slowly for these surface integrals to be non-zero. 

Consider initial data (say at t = 0) for a ffist order perturbation h^^^ that is of compact 
support (this will necessarily involve harmonics {l,m) corresponding to traveling waves). How will 
the second order perturbation sourced by this first order perturbation behave? Near infinity (at 
least at early times), h^"^^ must satisfy the source free linearized Einstein equation, i.e., the same 
equation as the ffi'st order perturbation. Hence the behaviour of h^"^^ near infinity should be the 
same as that of a ffist order perturbation. However, none of the first order normal modes decays 
sufficiently slowly to make a non- vanishing contribution to the surface integrals for the charges. 



e.g., a non-vanishing contribution to the energy would require r/ < 1 in (3.43) (with the lower sign 
choice) whereas we have seen that normal modes have r/ > 2.74. A non-vanishing contribution 
to the angular momentum requires rj < —1. Therefore /i^^^ does not behave like a normal mode 
at infinity. Furthermore, even the traveling waves decay too slowly to contribute to the surface 
integral for the angular momentum. So what linearized solution does /i^^^ behave like near infinity? 

Precisely the same issue arises for a Kerr black hole. Gravitational perturbations with I > 2 
decay too fast to contribute to the surface integrals for the energy or angular momentum. For Kerr, 
the resolution is that the Teukolsky or Hertz potential formalisms miss certain modes, specifically 
those modes that preserve the type D condition to first order. For Kerr, it has been shown that 
the only such perturbations correspond to deformations towards a nearby type D solution |31| . 
The nearby solutions are: the Kerr solution with different (M, J), the Kerr-NUT solution, and 
the spinning C-metric. The latter perturbations are excluded by asymptotic boundary conditions 
or regularity. Hence, for Kerr, one must add by hand the non-dynamical modes corresponding to 
infinitesimal variations in the mass and angular momentum of the black hole, which we can regard 
as / = and I = 1 perturbations respectively. Clearly these will decay at an appropriate rate to 
contribute to the surface integrals. 

This suggests that, in our case, the fall-off of /i^^) wih be the same as that of linearized modes 
that preserve the type D property. There are two classes of such modes: (i) modes that are locally 
gauge, i.e., locally of the form V(^r/^), and (ii) modes corresponding to a non-trivial deformation 
towards a type D solution continuously connected to NHEK. 
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Consider first the case that h^'^^ behaves asymptotically as a linearized mode that is locally 
gauge. By this we mean that, in a neighbourhood of the S'^ on which a boundary integral is 
computed, /i*-^^ is locally, but not globally, of the form V(^r7,y)|^ This is precisely what happens for 
Einstein gravity in AdSs, for example, where rj^ cannot be globally defined on the boundary. 
Could the same thing happen here? One might consider infinitesimal diffeomorphisms of the GHSS 
form e{(l))d/d(j) — e' {(j))rd / dr and, instead of taking e to be periodic in (j) (which would be globally 
defined), take e(0) = 0, which leads to a metric perturbation independent of t and (j). However, 
this has the effect of introducing a conical singularity into the metric near infinity (at the poles of 
the 5^), which does not seem appropriate. 

This "locally gauge" behaviour would arise from solutions that are obtained by identifications 
of the NHEK background (in the same way that the BTZ black hole is obtained as an identification 
of AdS^). Could one obtain a "NHEK black hole" by identifying the NHEK geometry in some 
way? Assuming any such identification acts only on the surfaces of constant 9, the possibilities 
have been well-studied [32] , and there appears to be no candidate free of pathologies such as conical 
singularities or closed timelike curves. 

Consider then, the second possibility, that h^"^^ behaves asymptotically as a linearized mode 
corresponding to a deformation towards a nearby type D solution. What solutions are there? 
Using Kinnersley's classification of type D solutions [33], the only such solutions appear to be: 
NHEK with a change in the angular momentum, the full (asymptotically flat) Kerr solution, or the 
near-horizon geometry of the extremal spinning C-metric. The latter has a conical singularity and 
so presumably must be excluded. 

It appears that the only candidate for a "Z = 1" mode, i.e., a mode contributing to the surface 
integral for angular momentum, is the perturbation that corresponds to a change in the angular 
momentum of the NHEK geometry ( J — > J -|- (5 J in ( |2.1[ )). This violates the GHSS fall-off condi- 
tions. Hence it would appear that, at second order, any perturbation with non- vanishing angular 
momentum is excluded by the fall-off conditions. 

What about the energy? One can attempt to obtain a solution with non-zero energy by taking 
a decoupling limit of the near-extremal Kerr solution at fixed temperature and angular momentum. 
An analogous decoupling limit of Reissner- Nordstrom was discussed in Ref. [31]- However, in the 
latter case, it was shown that, even with non-zero temperature, the decoupling geometry is simply 
AdS2 X S"^. We find that the same is true for Kerr: in Appendix [pj we show that the decoupling limit 
at fixed non-zero temperature leads back to the NHEK geometry. The explanation is presumably 
the same as in Ref. [M] , namely that the extreme Kerr black hole has a mass gap. 

It appears that the only regular modes with non-zero energy correspond to going to next order 
in the decoupling limit. This is probably equivalent (up to a SL(2,R) transformation) to retaining 
the next to leading order term in the near-horizon limit leading from extreme Kerr to NHEK. This 
clearly gives a solution fc^^^ of the linearized Einstein equation. However, it violates the GHSS 
fall-off conditions, indeed at the fully nonlinear level it amounts to considering an asymptotically 
flat black hole rather than its near- horizon limit. 

^"if it were globally a gauge transformation then it would give a vanishing boundary integral since, for a Killing 
field ^, Q^[h, g] is invariant under /i^,^ — > /ip,^ + V(^?;^) (even if 77 is a non-trivial element of the asymptotic symmetry 
group). 
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If correct, this implies that, if the first order perturbation has any non-zero energy (whether 
positive or negative) or angular momentum then at second order there will be a violation of the 
GHSS fall-off conditions This would be a satisfying conclusion: one does not have to worry 
about negative energy initial data, and the positive energy condition is redundant (at least in 
perturbation theory). What about initial data for a linearized gravitational field with vanishing 
energy and angular momentum? There certainly exists initial data with this property. We have 
seen that the normal modes have positive energy, so this data must involve traveling waves. With 
outgoing boundary conditions, the linearized theory predicts that these will disperse, leaving behind 
only normal modes, with positive energy. If this extends to the nonlinear theory then there still 
would be a problem since the final state would have to violate the fall-off conditions. It seems to us 
that the only solution is that, even though this initial data has vanishing energy, the two boundary 
integrals for the energy would be non-zero, but opposite in sign. Hence one would still obtain h^'^^ 
with the asymptotic behaviour just discussed, and thereby violate the fall-off conditions. 

This reasoning suggests that, at the nonlinear level, there are no non-trivial (i.e. non- isometric to 
NHEK) solutions of the Einstein equation that are continuously connected to NHEK, and satisfy the 
GHSS fall-off conditions (see Ref. [33] for a proof of a similar result for AdS2 x 5^). This may imply 
that the only solutions that satisfy the latter are related to NHEK by large gauge transformations. 
However, in this case, the dual CFT would consist purely of conformal descendents of the vacuum, 
which leads to a problem with modular invariance. Alternatively, there might be further solutions 
that are asymptotic to NHEK in the GHSS sense, but not continuously connected to it. If so, it 
would be interesting to find these solutions. 
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If one wanted to impose these fall-off conditions only at r ^ oo but didn't care what happened at r ^ — oo 
then one could always add to /i''^' an appropriate multiple of fc^^^ to arrange this, because k^^^ satisfies the linearized 
Einstein equation, and one is free to add to /i'^' (which satisfies the inhomogeneous equation ( |3.58[ )) any solution 
of the linearized Einstein equation. Note that the boundary integrals associated with any solution of the linearized 
Einstein equation, e.g. k^^\ must sum to zero. 
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Appendices 



A Shear-free null geodesies and NP tetrad for NHEK 
A.l NP quantities for NHEK in global coordinates 

NHEK is a Petrov type D geometry and its NP null tetrad is found by looking into the congruence 
of shear-free null geodesies |2l!. 

Geodesies are the paths that minimise the action associated with the Lagrangian 



1 dx^' dx" _ 5 



(A.l) 



where A is an affine parameter, and 5 = 0,1, respectively, for null and time-like geodesies. Since 
the NHEK geometry (2.1) is stationary, the energy E and angular momentum L of the particle, 

Pt = gtd + gtipif = E , p^ = gt^i + g^^ip = L , (A.2) 



are conserved in a geodesic motion, where = is the conjugated momentum. Equation (A.2) 
yields 



E-Lr 



" M2(l-hcos2 6l)(l-hr2) ' M2(l-hcos2 6i)(l-hr2) 4M2 sin2 6^ 

The Hamilton-Jacobi equation for the geodesic motion on a geometry g^y reads 



r{E — Lr) 



L 1 + cos2 



as 

9A 



1 dS dS 



Assuming a separation ansatz of the form 

1 



S=-5X + Et + L^ + Sr{r) + Se{e) 



(A.3) 



(A.4) 



(A.5) 



equation (A.4) for the NHEK boils down to 

2 



dSe 

de 



= m , 

V dr 



9 (l + COs2 6) 

Q{e) = A - l2 ^ . . . „ ^ - M^5 cos2 e 



n{r) 



n{r) = {E- rhf - (A + m2<5)(1 + r"^) ; 



(A.6) 



where A is the separation constant. Typically one has f'^ oc lZ{r) and 9"^ oc 0(0) and the conservation 
equations (A.3) give the remaining equations for t and 93. 

The shear- free principal null geodesies are found by requiring 9 (x Q{9) = for 5 = 0, which in 
our case requires 



A = L2 



;i + cos2 9oy 

4 sin2 00 



(A.7) 
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for a constant 9 = ^o- Moreover, these geodesies must also keep r, i.e., 

-.0 , o^ 9 (l + cos^ 6*0)^ 
n={E- rLf - (1 + r2)L2 ^ . . ^ 

4 sm t^o 

constant along the motion. Clearly this is possible for any r and 6q only if L = 0. The energy and 
angular momentum conservation equations (A. 3) then require 

E 



M^(l + cos^ OqY 



M^(l + cos^ eo)ip 



rE 



(A.9) 



1 + 1 + 

This Hamilton-Jacobi analysis then concludes that shear-free null geodesies have the tangent vectors 



1 



1 +r2 
1 



5t + 9, 



1 I 2 ' 
1 + 



(a, - (1 + r2)a, - r9^) 



(A.IO) 



4M2fi2(0) 

that we choose for the real vectors of the NP tetrad since they satisfy the appropriated relations 
in (A. 16). In particular, the normalisation factor for was chosen to satisfy the normalisation 
condition £ ■ n = 1. 



We can now check that (A.IO) are indeed null geodesic generators and in our way we find 
Carter's constant of motion for the NHEK. With the NP tetrad choice (2.3) we find the Weyl 
scalars (2.4) and NHEK is then Petrov type D. In such a spacetime, if we take = {t,r,6,ip) to 
be an affinely parametrised geodesic, k^^Vf^ku = 0, then |21] 



K = 2|^2r^/^(k-£)(k-n) -Q|k|2 

= 2|^'2r^/^(k-m)(k-m) - (Q- 1^' 



-^2/3 



is conserved along k if and only if a scalar Q exists which satisfies the equations 



DQ = D\^2\ 



-2/3 



AQ = A|^2 



'2/3 



SQ = 5*Q = 0. 



(A.ll) 



(A.12) 



In our case, from (2.4), one has |*I'2 



-2/3 



Af''/^ (1 + cos^ 6). The scalar Q = M^l^ satisfies ( |Al2| ) 
and we can then construct the two conserved Carter quantities (A.ll). This yields the pair of 
equations 

(1 + cos^ efe'^ = -L^ - [5 cos^ 6 - K) 

^ ' 4sin"'6' ^ 



(1 + cos^ efr'^ = {E- rLf - M'^{5 + K){1 + 



(A.13) 



whose RHS is, respectively, 0(0) and 7^(r) defined in (A. 6) if we identify A = M K. Carter's 
equations (A.13), combined with the energy and angular momentum conservation equations (A.3), 
reduce the finding of geodesies in NHEK to a quadrature problem. 

If we want shear-free null geodesies we demand 5 = and = which implies the relation 
(A.7). The radial equation then stays 



(1 + eos^ Ooff'^ = {E- rLf - (1 + r^)^^ 



(l + cos2 Oof 
4 sin^ On 



(A.14) 
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which can be independent of r only for L = 0. Inserting this condition in (A. 3) yields (A.9), and 
under the redefinition E ^ E (^1 + cos^ 9oj we finally confirm that shear-free null geodesies are 
those that satisfy 

E 



t=—^, r = ±E, 6 = 0, ^ = -T—^E. (A.15) 

These two geodesies give us the null NP vectors i and n as well as the Eddington-Finkelstein 
coordinates for the NHEK. The NP tetrad is completed with the introduction of the complex 
conjugate pair of vectors and m'^ as defined in (2.3). These are found requiring that the NP 

0, 

(A.16) 



tetrad satisfies the normalization and orthogonality conditions 

£m = £m = nm = nm 
££ = nn = mm = mm 
£n = l, mm=— 1. 

In terms of the NP tetrad, the metric components read 



0, 



(A.17) 



The 12 complex spin coefficients are introduced through linear combinations of the 24 Ricci 



rotation connection coefficients 7ca6 = ^(^c)^{h)'^ '^^{a.) 



r 



7311 
7243 



0, 
0, 



O" = 7313 
P = 7314 

isin^ 



0, 
0, 



7312 
= 7241 



V2M {I + cos^ ey 

i sin 6 



/2M(1 



7 



a 



0, A 

(7212 + 7342) = 



I cos I 



7242 
1 

2 
1 
2 
1 
2 



7244 



0, 
r 



1 



(7211 + 7341) = 0, 



(7214 + 7344) 
(7213 + 7343) 



2M2 (1 + cos2 9] 
cos(9-if2 



COS' 



e) 



sin 

cos 9 



(A.18) 



+ i cos 9) sin 9 

Their complex conjugates are obtained through the replacement 3 <-> 4 in jcab- From the Goldberg- 
Sachs theorem, K = cr = z^ = A = implies that the NHEK is Petrov type D (as it must be by 
construction). Moreover, e = implies that £ is affinely parametrised as it is indeed the case. 
The Weyl tensor 



C, 



R 



1 



puaf3 



+ ^{giiagui3 - g/iijgua) , (A.19) 



reduces to the Riemann tensor because (2.1) is Ricci flat. The 5 complex Weyl scalars ^'j in the 
NP formalism encode the information on the 10 independent components Cabcd of the Weyl tensor, 



^0 = 


— ^1313 = 




^1 = 


— ^^1213 = 




^2 = 


— C1342 = 




^3 = 


— Cl242 = 




^4 = 


— C2424 = 





(A.20) 
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For the NHEK these Weyl scalars are hsted in (2.4). 



The fundamental quantities in the NP formahsm needed to study perturbations are the spin 



coefficients hsted in (A. 18) and the directional derivative operators, 



D = , A = n^V^ , 5 = m^V^ , 5* = m^V^ . (A.21) 
A. 2 Master equation for NHEK in Poincare coordinates 

For completeness we write here the master equation for NHEK in Poincare coordinates. This is 



the counterpart of the global coordinate master equation (2.5). 



Let quantities with tildes denote Boyer-Lindquist coordinates of the full black hole solution 



(D.l), and take {r, y?} to be the Poincare coordinates describing NHEK. 

Bardeen and Horowitz define the near-horizon limit of the extreme Kerr solution by setting [1] 

r = a + \y, t=^, <^ = (^ + -^, (A.22) 
A ZaX 

where a is the extreme value for the Kerr rotation parameter, and taking the limit A ^ with the 
untilded quantities held fixed. The limit yields the near-horizon solution in Poincare coordinates. 
Taking this limit in the Kinnersley tetrad [13] , one finds that 

^ 2a2 a d Id 

y Ot dy y d^p ' 

1 1 f 2 d y'^ d y d 
A n ^ -^r- — a 



a'^{l + cos0)\ dr 2 dy 2 dip ) ' 
1 (9 , . (l + cos^g) d\ 

^/2a{l + icose)\de 2sm9 dp)' ^ ' ' 

Hence, by performing a boost before taking the limit, we can ensure that the tetrad remains well- 
defined and must therefore give a tetrad aligned with the principal null directions of the near-horizon 
geometry. 

Consider the Teukolsky equation for a field ^ of spin s in the extreme Kerr geometry [13] . Let 

^ = /(r, y, GM = ! (a*, ^ ' ^' - ^) ' (A-24) 
Plugging this into the Teukolsky master equation [13 , and taking A ^ 0, we find that it becomes 



dlf -^drd^f+(2-- sin2 9 - d^f - y~^^dy [y^'+^dyf) - -de (sin Odef) 

y \ 4 sm^ 6 J ^ ^ ^ sm 6' 

-is C^^^ + cos e] d^f - —drf + (s^ cot^ e-s)f = G (A.25) 



4a^ 

y2 

—is I 

sin^ 8 J y 

This is the master equation governing perturbations of the near-horizon Kerr geometry written in 
Poincare coordinates. 

We separate variables by setting 

f{T,y,e,p)=F{T,y)S{e)e'^^. (A.26) 
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Equation (A. 25) separates into an angular equation given by (2.7) with C = m/2 and into 
4a^(s + im) 



y 



drF - y-^'dy {y^'+^dyF) + A, 



(A.27) 



where A^^ is the constant of separation already discussed after (2.7). 



B Phase and group velocities 



In this Appendix we give some details of the analysis done in subsection 2.5.2 

To discuss travelling waves and their phase and group velocities we need the next-to-leading 



contribution to the asymptotic behavior (2.27). This is obtained applying to (2.26) the transfor- 
mation law z ^ 1/(1 — z) of the hypergeometric function and its asymptotic expansion for large 
radial distances yielding, 



(s) 
Imuj 



2(l+»7)/2p(^_ ^)(-v±g±i7r(/3-a-a)/2g-i±2ln|rh-^_^^^ 



2<;c^ 1 



_ 1-2, 



I I 1— r/ r 



(B.l) 



where the amplitudes and are defined in (2.28). 
Introduce the quantities 



S±{r) 



k±{r) 



exp 



I ±-?7ln r + —^32 

2 1 + T]'^ r 



7] = ir] , 



drS±{r) 
S±{r) 



(B.2) 



Here, S±{r) encodes the radial contribution to the wave propagation and, in a WKB approximation 
to the traveling waves, k±{r) is the effective wavenumber of the wave. The superscript ('— ') is 
in correspondence with the amplitude (D^). The phase velocity of the traveling wave is then 



UJ 



2ujr 



k±{r) 



while the group velocity is 



s I duj 



1 {l + rf ) 

2 mrj^ s 



(B.3) 



(B.4) 



For 2 < / < 20 and the values of \m\ < I that yield 77 > we have checked that, in the denominator 
of w^, mry =p s is positive for m > and negative for m < 0. 

Both at r = ±00, depending on whether we choose the or the contributions in ( |B.1[ ), 
we can have the combinations for the sign of the phase and group velocities displayed in Table |3] 
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C Hertz map between Weyl scalar and metric perturbations 



In this Appendix we derive the map hf^^{^) between the Weyl scalars ^0,4 (that satisfy the decou- 
pled Teukolsky equations) and the metric perturbations h^^, as well as the map Afj,{(j)) between 
the NP complex scalars perturbations ^0,2 and the Maxwell perturbed vector potential A^. We 
follow |18] and recover the results of [151 Ull [13 US ES] but we take the opportunity to emphasize 
the importance of these works (not so well-known in the community) and to give some details not 
presented in the original articles and to also pinpoint small typos in some of these works that have 
propagated in the literature. We start by reviewing Wald's work |T5] in the next subsection. Then 
we apply it to get hfj_^{"^) (subsection |C.2 ), and to find ^/i(</>) (subsection C.3). A similar analysis 



could be done to obtain the map for the Weyl fermionic perturbations. 
C.l Problem statement. The Hertz potential map 

Let us start by briefly reviewing the seminal work of Teukolsky |12| [13]. Suppose we wish to 
solve the perturbation equation e{h) = where e is a linear differential operator and h is the field 
perturbation on which e acts. For example, e can be the Maxwell operator describing electromag- 
netic perturbations A^, [AE(^a)]^ = 0, i.e., V ^F^^ = 0; or it can be the gravitational operator 



describing gravitational perturbations hap, [AG(/ia/3)]^^ = written in (3.42 1 . 
Suppose now that: 

• a new variable ^ = L^^{h)^ where is a linear differential operator, has been introduced 
(these are e.g., the NP complex electromagnetic scalars </>o,i,2 or the Weyl scalars ^'o,---,4); 

• a linear partial differential operator D has been found such that for all h one has 

I])e{h) = OL^{h) = 0^ , (C.l) 
where O is another partial differential operator. 

Then 

e(/i) = O(^) = 0. (C.2) 







c- 


D- 


c+ 


D+ 
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> 
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m > 
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> 




m < 


> 


< 
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< 



Table 3: Phase and group velocities for the possible amplitudes choices in the asymptotic solution 
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This is the main result of [I2| I13j. who found the variable ^ and operators L\j,, D and O, as well as 



the well-known decoupled Teukolsky equations ( C . 1 ) , that describe the problem of electromagnetic, 



Weyl fermionic and gravitational perturbations in the Kerr black hole. 

So assume that we have carried on the previous steps that lead to (C.2) and that furthermore 
we solved it and have a solution for ^. The next question is how to get the original perturbation 
h from the knowledge of the scalar perturbation ^, i.e., the unique map h = This issue has 

been addressed by Cohen and Kegeles [151 llZj i Chrzanowski [16] and later Wald |T8] proved 
rigorously and in a few lines their results. In the sequel we review Wald's proof |18|. 

Start by recalling the notion of adjoint of an operator. Let O be a linear differential operator 
taking a scalar, vectorial or tensorial field into another similar field. Then there is an unique adjoint 
operator such that 

(O^) - (o^^h) ^ = V^s^ (C.3) 



for arbitrary fields ^ and ^'h, and where V^s^ is a total divergence. 
Wald's theorem states the following. Assume that: 



• the identity (C.l) is satisfied for the linear differential operators D, e, O, L^; 

• ^'h satisfies O^^^h = 0, where ^'h is called the Hertz potential; 
Then 

• Dt^'H satisfies (D^'^'h) = 0. 

• Moreover, this in particular also implies that if e is self-adoint, = e, then 

h = D^H is a solution of e{h) = . (C.4) 



and provides the map we are looking for. 

The proof of this result is short and simple [18]. Taking the adjoint of Oe = OLx^, one has 
^tpt = l\,0^ . Applying these operators to ^'h one gets e^'D'^'^'H = after using the assumption 
Ot^'H = 0. Moreover, if et = e, it trivially follows that e (O'^^'h) = 0, i.e., h = Dt^fH is a solution 
of the initial perturbation equation e(/i) = 0. 

Finally note that the Hertz potential construction reviewed here applies to geometries where 
no energy-momentum tensor is present. For a discussion of the method when this is not the case 
as well as for the second order perturbation analysis we ask the reader to see |25| . 



C.2 Application: the map h^u{^) for gravitational perturbations 

Our starting point are the decoupled Teukolsky equations for the perturbed Weyl scalars ^'g^^ and 
'^^^\ namely equations (2.12)-(2.15) of [13^, which are written in (C.16) for s = 2 and (C.17) for 
s = —2. These can be written as (to make the connection with the nomenclature of the previous 
subsection straightforward) 



(C.5) 
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where we use the superscript (1) to denote a perturbed quantity (otherwise it refers to an unper- 
turbed quantity), and 



Ogo = - 3e + e - 4/) - p)(A + ;u - 47) - (5 + 7f - a - 3/3 - 4r)((5 + vr - 4a) - 3^'2 , 
Og4 = (A + 37 - 7 + 4;u + iu)(L' + 4e - p) - - T + ;5 + 3q + 47r)(5 - T + 4/3) - 3^2 : 

and 

(5 + 7f - a - 3/3 - 4r) \{D - 2e - 2p)e^V^ - (5 + vf - 2a - 2/3)£^r 



(C.6) 



+{D - 3e + e - 4p - p) [((5 + 27f - 2(5)6^'m''^ - {D -2e + 2e-p) 



(C.7) 



(A + 37 - 7 + 4/i + p) ((5 - 2r + 2a)n(^m'') - (A + 27 - 27 + ji)m^'rff 

+(^ - r + /3 + 3a + 47r) [(A + 27 + 2/i)n(^m'") -(5-T + 2'p + 2a)n^n'' 

Following Wald's procedure, the Hertz potential \I'h is introduced to be such that it satisfies 
the equation Oq^'h = 0, that is 

[(A + 37 - 7 + Ji){D + 4e + 3p) - (J + ;5 + 3a - t)(5 + 4/? + 3t) - 3^2] ^Hq = , 
[(L> - 3e + e - p)(A - 47 - 3p) - (5 - 3/3 - a + 7f)(^ - 4a - 37r) - 3^'2] ^'h4 = , 

where, to obtain the adjoint of (C.6), we used the relations^ 

At = - 



(C.8) 



(jt = _(5 + /3_a-r + 7f), 



ft 



_(A - 7 - 7 + p + /i) 
(5 + /3-a-r + 7r), 



(C.9) 



and the well known property {AB)'^ = B^A^. Equations (C.8) can be written, respectively, as 
dClsl ) with s = -2 and (jO^ with s = 2. 

Since the gravitational perturbation operator e = [AG{hap)]^j^ is self-adjoint, Wald's theorem 
tell us that the map (C.4| between the Hertz potential ^'h and the metric perturbations h^iy is 
given by /i^j. 



2Re 



I.e. 



h 



IRG 



c.c. 



i,ORG 



|£(/,m^) [{D + 3e + e - p + p){6 + 4P + 3t) + {6 + 3P - a - T - Tf){D + 4e + 3/))] 

-l^,^y{5 + 3/3 + a - t)(5 + 4/3 + 3t) - m/,m^(D + 3e - e - p)(D + 4e + 3p)|*Ho 
^ni^fjfiy) [(5 + /3 - 3a r 7r)(A - 47 - 3/i) (A - 37 - 7 p - Jl)(5 - 4a - 37r)] 
-n^nu(5 - /? - 3a -I- 7r)((5 - 4a - 37r) - m^my{/S. - 37 -|- 7 -|- /i)(A - 47 - 3/x) >\I'h4 + c.c. , 



(C.IO) 



^^To get (C.9 1, introduce the internal product (^'HjO'i/') = d^x-^—g '^YiOtji. For the directional derivative 
operators O = fij-^^V^j it then follows, after integrations by parts and use of V^^/^ — 0, that {^'HjC'f/') = 



-(t/^.O^^h) + £" VpS" with 



and 'Vf^s'^ being the short notation for a total divergence 



contribution. To compute V^e^^j use the well-known relation between covariant derivative of the NP tetrad and the 



spin coefficients [5T], V^jef^)^ = e j^j^ate and relations ( |A.18[ ) 
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where to get the adjoint of (C.7) we used again (C.9), and c.c. stands for complex conjugate. The 



first of these relations gives the metric perturbations in the ingoing radiation gauge (IRG), while 



the second provides the map in the outgoing radiation gauge (ORG); see (3.39). The first relation 



agrees with the results of |15| [T6l IT71 118]. and is equation (3.41 ) in the main body of the text. The 
outgoing radiation map corrects typos in the relation of Table 1 of [16] that have propagated in 
the literature. 



C.3 Application: the map A^{(f)) for electromagnetic perturbations 

We begin our discussion with the decoupled Teukolsky equations for the perturbed electromagnetic 



NP scalars <^q^^ and namely equations (3.5)-(3.8) of [13], which are written in (C.16) for s = 1 
and (C.17) for s = —1. These can be written as 



/(I) 



where, again, we use the superscript (1) to denote a perturbed quantity, and 

OEo = {D-e + e-2p- p){A + - - {6 - ^ - a - 2t + 7f)(6 + n - 2a) , 
= (A + 7 - 7 + 2/i + JI){D + 2e - p) - (6 + a + P + 27T -t){6 - T + 2p) , 



(C.ll) 



(C.12) 



and 



{6-a- p + Tf- 2T)t - {D -e + e-2p- p)mi' , 
D{^^ = (A + 7 - 7 + 2^ + - (6 + a + ]3 + 2Tr - T)n>' . 

The Hertz potential ^'h is introduced to be such that it satisfies the equation Oe^h 



IS 



[(A + 7 - 7 + 7x)(L» + 2e + p) - ((5 + /? + a - t)((^ + 2/3 + r)] ^Ho = , 
[{D-e + e--p){A-2-i-p)-{5-f3-a + -2a- ^)] = , 

where, to obtain the adjoint of (C.12), we used (C.9) and {AB)'^ = B'^ A'^ . Equations 



(C.13) 
0, that 

(C.14) 



C.14 



are, 



respectively, equations (C.18) for s = —1 and (C.19) for s = 1. 



Since the Maxwell perturbation operator e = [AeC^q)]^ is self-adjoint, Wald's theorem tell us 



that the map (C.4) between the Hertz potential ^'h and the vector potential perturbations is 



given by = 2Re 



I.e., 



A^RG ^ + 2/3 + r) + m^{D + 2e + p)\ ^h,, + cc. 

aORG 



[-m^(A - 27 - /i) + n^,{5 -2a- vr)] ^'Ha + cc. , 



(C.15) 



where to get the adjoint of ( C.13 ) we used again ( C.9 ) . The first of these relations gives the Maxwell 
perturbations in the ingoing radiation gauge, i^A^^ = 0, while the second provides the map in the 
outgoing radiation gauge, n'^j4^ = 0. The first relation agrees with the results of [151 Ull 113 118] . 
The outgoing radiation map agrees with [1^ and corrects typos in the relation of Table 1 of [THI . 
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C.4 Hertz potential for spin s in a Ricci flat Petrov type D geometry 



For an arbitrary Petrov type D solution we can write the decoupled Teukolsky equations in a 
single equation that depends on the spin of the field |13] , and the same happens for the decoupled 
equations for the Hertz potential if the the solution is furthermore Ricci flat. 

For positive spin and negative spin, Teukolsky's decoupled equations are, respectively, [13] 



[D - (2s - l)e + e - 2s/9 - ^] (A + /i - 2s7) - [<5 + vf - a - (2s - 1)/? - 2sr] {5 + tt - 2sa) 



-3s («- 2 1(^-1)^2 



^W=47rT(,), for s = +2, +1, , (C.16) 



[A - (2s + l)j - ^ - 2s^i + Jl] {D - 2s€ - p) - [6 - T + P - (2s + 1)q - 2s7r] {6 - t - 2s/?) 



+3s ( s + - ) {S + 1)^2 



= 4vrT(^-,, for 



-2,-1, --,0. (C.17) 



The relation between the nomenclature used here and the original notation of Teukolsky [13^ was 
already displayed in Table [T] and associated discussion. 

The Hertz potential in the ingoing radiation gauge describes negative spin s = — 2,— 1,— 2 

(s) 1 

field perturbations, while in the outgoing radiation gauge ^'^ describes positive spin s = +2, +1, + 2 
perturbations. They are, respectively, the solutions of the scalar equations 



[A - (2s + 1)7 - 7 + /I] - 2se - {2s + l)p]- [6 + P - (2s + l)a -f][S- 2s/3 - (2s + l)r] 



^\f = 0, for s = -2,-l,--,0, 



(C.18) 



[D - (2s - l)e + e - p] [A - 2s7 - (2s - l)p] -[5- (2s - l)p - a + Tf] [6 - 2sa - (2s - l)7r] 



-3s( s - - ) (s-l)^2 



W = 0, for s = +2,+l,+-,0, 



(C.19) 



The conjugate Teukolsky perturbations ^'^^^ to these Hertz potentials ^'^ ^ can be read from Tables 
□ and 12 



For the NHEK geometry written in global coordinates the NP tetrad is written in (2.3 ), the spin 



coefficients are listed in (A. 18), and the directional derivative operators can be read from (A. 21). 



Using this information in (C.16)-(C.17), and (C.18)-(C.19) we get the master equation (2.5). 



D Decoupling limit of near-extreme Kerr. Mass changing modes 

In this appendix we show that a decoupling limit of the near-extremal Kerr black hole yields the 
NHEK geometry. We follow [31] where decoupling limits like the one we take were first discussed 
for charged non-rotating solutions. 
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The Kerr black hole solution in Boyer-Lindquist form reads, 

ds^ = ^ — d? -^dr^- ^de^ 

(f2 + a2)2 _ Aa2sin2 6' A 

(f^ + a^Y -Aa^sin^e f ~ a (a^ + - A) V , ^ 



(f2 + a2)2 - Aa2sin2 



with 



A = (f — ) (f — f+) , E = + cos^ Q , 

f± = l {fpAE + Ip^ epAE-^ + A.J^ ± \ji%AE^f^AE'^ + AJ . (D.2) 



We used the fact that in four dimensions the Planck length £p is related to Newton's constant G 

-P 



by = G, and we defined the excitation energy above extremality as 



AS=!?1^, with a = 4-, (D.3) 
Ip m 



where J and M = m/i^p are the ADM angular momentum and mass of the Kerr black hole. 
The black hole temperature is 

Th = ' 

47rr-(_ (^r^ + a^j 

from which follows that near-extremality the relation between the excitation energy and the tem- 
perature is 

~ SvrVtT^^P. (D.5) 

Typically, the energy of a quantum of Hawking emission is of order Th- When this energy is of 
order or greater than the available energy above extremality, Th > AE, the semiclassical analysis 
of the black hole thermodynamics breaks down. This occurs at an excitation energy of order 

E,^, - -A- • (D.6) 

We now want to take a decoupling limit where 

£p ^ , with (Tff , J) fixed. (D.7) 

In this limit the excitation energy AE vanishes and the gap energy -Egap goes to infinity. 

In the decoupling limit, and after introducing the new radial and azimuthal coordinates {Ujip): 

f = f+ + 2£lU, 4> = ij + ^, (D.8) 



the Kerr geometry (D.l) reduces to 
ds^ 



ip 



J2 



d? + jj^^yf^^ + dO^ + (.) {d^ + {2.TH + ^) , 

(D.9) 
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with 02((9) and h?{e) defined in pl2| ). 

Finally if we introduce the new time, radial and azimuthal coordinates (r, y, ip) through the 
transformations, 



1 



U 



27tTh 
_ ttTh J 

y 



arctan (t — y) + arctan (r + ; 

(y-i)'-T 



ip = (p + arctan 



1 



y 



+ arctan 



T + 1 

y 



the decoupling geometry (D.9) reduces to 



4 



(D.IO) 



(D.ll) 



We recognize this geometry as the NHEK solution written in Poincare coordinates. A final trans- 
formation between the Poincare coordinates (r, y, 9, if) and the global coordinates {t, r, 6,(1)) |lj. 



sm 



t Vl + r"^ + r 



cos 



cos 



2 j 



+ ln 



cos i + r sin t 



1 + sin t \/l + r'^ 



(D.12) 



takes (D.ll) into the NHEK (2.1) written in global coordinates. Therefore, at the classical level, 
when we take the decoupling limit (D.7) of the near-extremal Kerr geometry, one gets a geometry 
that is independent of the temperature Th and that is precisely the NHEK geometry. 

An important observation is that going to the next-to-leading order in the ip expansion, a path 
similar to (D.7)-(D.12) yields the next-to-leading order contribution to (D.ll) or (2.1), that we call 
hfj_^. This hfj_i, can be considered as a perturbation to the NHEK since it satisfies the linearized 
Einstein equations (with h ^ 0) hy construction. Asymptotically, h^^ goes as a power of higher 
(in all components) than the GHSS boundary conditions [2j- 

To interpret physically these perturbations we compute the relevant conserved charges asso- 
ciated to them as defined in (4.62). One finds that the energy of these perturbations is finite, 
Qdt[h,g] 7^ 0, while the U{1) charge vanishes, Qd<f>[h,g] = 0. These perturbations therefore in- 
crease the energy of the solution while leaving its angular momentum unchanged. They correspond 
thus to perturbations that take NHEK away from the extremality state. Since the full series expan- 
sion in reconstructs near-extreme Kerr, these perturbations actually take the NHEK geometry 
into near-extreme Kerr black hole. 
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